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$H ' The electronic states near a surface or a domain wall in the p-wave superconductor are studied 

CL|' 

r- { ■ for the order parameter of the form p^^ztipy-wave, which is a unitary odd-parity state with broken 

"^ . time-reversal symmetry. This state has been recently suggested as the superconducting state of 

C^ , Sr2Ru04. The spatial variation of the order parameter and vector potential is determined self- 

consistently within the quasi-classical approxinration. The local density of states at the surface 

i-rt \ is constant and does not show any peak-like or gap-like structure within the superconducting 

Ch ' energy gap, in contrast to the case of the d-wave superconductors. The influence of an external 

fj ' magnetic field is mainly observable in the energy range above the bulk gap. On the other hand, 

there is a small energy gap in the local density of states at the domain wall between domains 

OQ ■ of the two degenerate Px+'ipy-'^au\e and p^,— ip^-wave states. 

> 
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Sr2Ru04 is the first superconductor with layered perovskite structure, which does not contain 

copper .Id) Although the structure is identical to that of some of the high-temperature superconduc- 

^ ' tors, the transition temperature is rather low, rc=1.5K. There is a clear difference in the electronic 

structure, since Sr2Ru04 is a good metal and even a Fermi liquid in its stoichiometric composition. 

K^ , Band structure calculations in good agreement with the de Haas-van Alphen measurements show 

j^ ■ that this compound has three Fermi surfaces originating from the three 4(i-t2g-orbitals of Ru^'^.BQ'q) 

There is growing experimental evidence that the superconducting state is unconventional (non-s- 

wave). Examples are the absence of a Hebel-Slichter peak in 1/TiT of NQR-measurementa3 and 

the sensitivity of Tq on non- magnetic impurities .Q' 

It was suggested that the superconducting state has odd-parity (spin triplet) pairing.QEHlJ'&EiP 
There is a certain similarity with ^He considering the correlation effects (superfluid ^He has p- 
wave pairing) .LP Furthermore, there is a series of related compounds such as SrRuOs which are 



ferromagnetic suggesting that ferromagnetic spin fluctuations are probably enhanced in Sr2Ru04 
and mediate odd-parity, spin triplet pairing£3l3'E3' The recent discovery of intrinsic magnetism in 
the superconducting phase by /uSR experiments indicates a pairing state with broken time reversal 
symmetry.EJ) Symmetry considerations lead to the conclusion that this would only be possible for 
an odd-parity state.llB' A very strong support for odd-parity pairing comes also from the Knight 
shift data in the ^^0-NMR measurements which demonstrate the absence of any reduction of the 
spin susceptibility in the superconducting statcllj^ The superconducting state compatible with all 
of these experiments is given by d{k)=z{kxiLiky) S3 

The presence of three electron bands forming the Fermi liquid state leads to the question whether 
all of them contribute to superconductivity on equal terms. Symmetry considerations show that 
it is possible to separate the orbitals {4:dyz,4:dzx} (forming two Fermi surfaces) and 4:dxy (forming 
one Fermi surface). Thus it is theoretically possible that superconductivity predominantly appears 
in one of the two sets of orbitals, while the other participates passively only through induced su- 
per conductivity.113) There are several supporting experiments for the scenario of "orbital dependent 
superconductivity". The specific healllM' and NQR experiments^ indicate a considerable residual 
density of states in the superconducting state which could be attributed to the passive electron 
bands. Furthermore, the analysis of the London penetration depth extrapolated to the zero tem- 
perature is compatible with taking only one of the two subsets of orbitals into account for the 
superconducting state.L3 The present experimental situation is compatible with a superconducting 
state dominantly in the 4(i2:j^-orbital. The recent prediction of a square form of the vortex lattice 
aligned with the crystal lattice axis was confirmed experimentally and is consistently explained 
with the assumption of superconductivity in the 4(ixi/-orbital as well.o'ES' It seems therefore to 
be justified to approximate the electronic structure by a single band model consisting a single 
practically cylinder-shaped Fermi surface. 

The theoretical study of boundary effects in high-temperature superconductors with d-wave 
Cooper pairing have turned out to be very fruitful for subsequent experimental investiga- 
tions. The bound states cause unusual features in the Andreev reflection and tunneling spec- 
troscopy£iHH3E2IB3'E3'E3'E3) Recently characteristic tunneling conductance properties were examined 
theoretically for the triplet pairing states.E§L2P Conductance peak features related to the bound 
states are very sensitive to the angle of the incidence in this case, in contrast to the d-wave super- 
conductor. In this paper we investigate the electronic states near the surface and compare them 
with those found near domain walls. Domain walls appear due to the fact that the superconducting 
state is two- fold degenerate {px+ipy- and px—'ipy-wave state) such that two types of domains are 
possible. Because of the discrete degeneracy the separating domain walls are very well localized and 
constitute a region of spatial deformation of the order parameter giving rise to subgap quasiparticle 
states analogous to the surface. 



§2. The Quasiclassical Formulation 

In this paper we study boundary effects in the p-wave superconductor, determining the order 
parameters self-consistently. For this purpose we use the quasi-classical Green function formalism 
developed by Schopohl, which, for example, has been applied to the vortex problem.L3'EiP We adapt 
this scheme to the boundary problem and apply it to both the case of the surface and the domain 
wall. 

For simplicity we assume that the superconductor is two-dimensional with a single cylindrical 
Fermi surface which is a reasonable first approximation to the band of the 4da;y-orbital, and that 
the superconducting order parameter has a p-wave symmetry. The gap matrix for the p-wave state 
is given by the d-vector defined as 

A{r,r') = idir,r')-aay, (2.1) 

where ai{i=x,y,z) is the Pauli spin matrix. We restrict ourselves to the case of 

d(r,r')=(0, 0, A(r,r')) type, so that the Bogoliubov-de Gennes equation reduces to a 2 x 2 matrix 

form, 

^^/5{r-r')ho{r') A(r,r') \ /n,(r')\ ^ / uiir)\ 

Here A(r,r') is the p-wave pair potential, Ei is the l-th energy eigenvalue, and 
/io(t')=[— iV+eA(r)]^/2m— i^F is the kinetic energy of an electron measured from the Fermi energy, 
with e (e>0) and A as the charge of the electron and vector potential, respectively. For convenience 
h is taken to be unity throughout this paper. 
Introducing the Andreev approximation (i.e. separating the rapid oscillations) as 

fu^ir)\f^iik,,r)\^,^,^^ (2.3) 

\vi{r) J \vi{kF,r) J 



we obtain the Andreev equation, 

/ -ivF ■ V + evp ■ A{r) A{kF,r) \ f ui{kF,r)\ _ / ui{kF,r)\ 

V -A*{-kF,r) WF-V + evF-A{r)J \vi{kF,r) J ^ \vi{kF,r) J ' 

where fcp and vf are the Fermi wave number and the Fermi velocity, respectively. The fcp de- 
pendence of A{kF,r) represents the symmetry of the order parameter and r is the center position 
of the Cooper pair. Due to the p-wave symmetry (i.e. A{—kF,r)=—A(kF,r)), eq. ( p. 4] ) has the 
same form as the Andreev equation in the case of singlet pairing. The corresponding Eilenberger 
equation is given by 

' iujrn - evF ■ A{r) -A{kF,r) 



it>F ■"^g{kF,uJm,r) 



A*{kF,r) -iLOm + evF ■ A{r) J 



(2.5) 



where g{kF,uJm,'r') is the quasi-classical Green function in a 2 x 2 matrix form and uJm=T^T {2m+l) 
{m: integer) is the fermion Matsubara frequency.Ej'E3'E3) For simplicity the Boltzmann constant is 



g{kF,UJrn,r) 



taken to be unity. The equations for the matrix elements of the quasi-classical Green function are 
expressed as 

g{kF,uJm,r) \f{k-p,ujm,r) \ 

-if{kF,uJm,r) -g{kF,uJm,r) J ' _ 
vp ■ Vg{kF,ujm,r) = A*{kF,r)f{kF,u;ra,r) - A{kF,r)f{kF,uJm,r), 

{um + ievF -Air) + -vp ■V)f{kF,uj.^,r) = A{kF,r)g{kF,uJm,r), 

{uJm + ievF -Air) - -vf ■V)J{kF,ujm,r) = A*{kF,r)g{kF,uJm,r). (2.6) 

These equations contain the following symmetries, 

g{-kF,uJm,r) = g*{kF,uJm,r), g{kF,-u;m,r) = -g*{kF,u;m,r), 

f{-kF,uJm,r) = -J*{kF,ujjn,r), f{kF,-ujm,r) = f{kF,ujm,r). (2.7) 

Using the following transformation, we can solve the Eilenberger equation in a simple wayJ3 

1 - a{kF, Wm, r)b{kF, LOm, r) 



g{kF,iOm,r) 

f{kF,UJm,r) 
J{kF,UJm,r) 



1 + a(fcF, Wm, r)b{kF, ujm, r) '' 
2a{kF,uJm,r) 

1 + a{kF,uJm, r)b{kF, lOm, r) ' 
2b{kF,uJm,r) 

1 + a(fcF, Wm, r)b{kF, lOm, r) '' 



(2.^ 



where o and b satisfy 



Vf ■ Va{kF,uJra,r) = A{kF,r) - A*{kF,r)a^{kF,uJm,r) -2 w^ + ief f -Air) a{kF,uJm,r), 
Vf ■Vb{kF,uJm,r) = -A*{kF,r) + A{kF,r)b'^{kF,(jJm,r) + 2 cOm + ieVF ■ A{r) b{kF,(jJm,r). 

(2.9) 

Let us now consider the case of a surface perpendicular to the x-direction. We assume that quasi- 
particles are specularly reflected at the surface. We solve eq. (2.9) along the classical trajectories as 
shown in Fig. 1(a) where the quasiparticle moves with the momentum close to A;fi from A to B and 
with momentum /uF2 from B to C. Since the surface is specular and translational invariant along 
the y- and z-direction, the incident momentum along the surface is conserved. We then match the 
two solutions using the boundary condition of the quasi-classical Green function at point B given 
byi'§) 

g{kFi,LOm,B) = g{kF2,u;m,'B), (2.10) 

which for a and b means. 



a{kFl,UJm,B) = a(fcF2,Wm,B), 6(fcFl,Wm,B) = 6(fcF2,Wm,B) 



(2.11) 



Since the system is translationally invariant perpendicular to the x-axis, the quasi-classical Green 
function or a{kF,uJm,r) and b{kF,uJm,f) depend only on x. Thus eq. (2.9) can be rewritten as 



-—a{kF,iOm,x) = < A{kF,x) - A*{kF,x)a^{kF,iOm,x) - 2 

ax VFx 



LVm + ievFyAy{x) a{kF,uim,x)>, 



-—b{kF,UJm,x) = < -A*{kF,x) + A{kF,x)b'^{kF,UJm,x) +2 Um + ieVFyAy{x) b{kF,iOm,x) \. 

dx VFx[ L '^ '^ J J 

(2.12) 

where vfx and VFy are the x and y component of the Fermi velocity, respectively. The initial and 
boundary conditions for eq. (2.12) are given by 

A(fcFi,oo) 



a{kFl,UJm,X = oo) 



'uJm + ieVFyAyioo)]"^ + |A(feFi,Oo)p + UJm + ievFyAy{oo) 
A*(fcF2,00) 



b{kF2,i^m,X = oo) = , 

^[wm + ievFy Ay {oo)]"^ + \A{kF2, oo)p + uJm + ievFyAy{oo) 

a{kFi,UJm,X = 0) = a{kF2,iOm,X = 0), 

b{kFi,UJm, X = 0) = b{kF2,^m, X = 0). (2.13) 

The position-dependent order parameter can be determined by the quasi-classical Green function. 
The gap equations for the px- and py-wave (i.e. A{kF,x)=Ax{x)cos9i:+Ay{x)sm6k) are given 
(using eq. (2.7)) aM'^ 
(Ax{x)\ w- /2F,cosefc\ 



1 



2 COS 9i 



■KTVp V - dOk 

G<ijJm<!jJC 



/{Ok, UJm, x) + f*{9k,UJm,x] 



2cOSek\f{9k,UJm,x)+ f {ek,u;m,x) - /{tT - 9k,UJm,x) - f {tT -ek,UJm,x) 



2sin( 



f{9k,UJm,x) + f*{9k,L0m,x) + f {lT - 9k, LOm, -x) + f* {iT - 9^, Um, -x) 



K, 



iog^+ E 



(2.14) 



Tc ^ /^ m-1/2 

0<m<wc/27rT 



where 9^ is the angle of the Fermi momentum measured from the k^-axis and ujq is the cutoff 
energy. We write A^; and Ay for the order parameter components corresponding to the px- and 
Py-wave, respectively. Generally, they are complex numbers. Tq is the superconducting transition 
temperature which we assume is the same for both components. 

The supercurrent due to the broken time-reversal symmetry is a general property, which was 
previously studied in the framework of Ginzburg-Landau (GL) theorySEH* Here we investigate it 
within the quasi-classical approximation. The current density along the y direction can be expressed 
by the quasi-classical Green function, 



1 f^ 
Jy{x) = -evFN{0)T E — / d9 k sin 9 k{-i7r)g (9 k, UJm,'- 



|<^m|<<^C 



0<Ulm<<^C 



2evFN{0)T Y^ / d9ksm9klm. g{9k,uJm,x)+g{TT -9k,uJm,x) 



(2.15) 



where the symmetry relations of eq. (2.7) have been used.cS'cJ N{0) is the normal state density of 
states per unit volume at E-p and n is the density of electrons, whereby both A^(0) and n include 
up and down spin electrons. 

The magnetic field and vector potential are calculated using the Maxwell equation. The magnetic 
field and the vector potential can be determined as follows: 

rx rx rco 

B^{x) = -/i / dxJy{x), Ay{x) = /i / dx'B^{x') - /i / dx'B^{x'), (2.16) 

Jo Jo Jo 

with /u as the permeability. Equation (2.16) is chosen to satisfy the boundary conditions of Bz{0)=0, 

Ay{oo)=0. Now we can solve the gap equation (2.14) iteratively by using eqs. (2.8), (2.12), (2.13), 

(2.15) and (2.16) until the self-consistency is achieved. 

§3. Quasiparticle Properties at the Surface 

3.1 Self- consistent solution 

Before examining the boundary effect, let us study first the bulk case, where / and / in the bulk 
are given by 

/(6'fc,u;„^,bulk) = 7*(6'fc,u;„,bulk) = Z' ' ^ , 

y^UJm^ + |A(6'fc,bulk)|^ 

A(6lfc, bulk) = A^(bulk) cos Ok + iAj^(bulk) sin 9k. (3.1) 

Here we have assumed real values for A^. and Ay, since we are interested in the px+ipy-wave state. 
Substituting eq. (3.1) into the gap equation, we can obtain the bulk solution at T=0: 

A^(bulk) = AJbulk) = A(0) = 2a;ce ^^ , Tq = ^e ^p , — ^ = vre'T ~ 1.76. (3.2) 

vr Tc 

Here 7 is the Euler's constant: 7=0.57721- • •. This solution is similar to the s-wave case, since 
[A(0fc,bulk)| is a constant on the Fermi surface. If one of the two order parameter components is 
absent, the solution becomes 

A^(bulk) = 0, Aj,(bulk) = A,i„gie(0) = 4u;ce"^i+^\ 

2e'^ujc --J- Asingic(O) 1 _^ , . 

Tq = e ^p , = 2e 2 7re ^ ~ 2.14. (3.3) 

vr Tc 

Note that Asingie(O) is larger than A(0). This is due to the fact that there are gapless points on the 
Fermi surface for the single existing solution and that the denomitators of / and / become small. 
Though A(0) is small, its solution is more stable, since it opens a gap everywhere on the Fermi 
surface JI3) 

We now turn to the surface problem. As in the case of the d-wave pairing, boundary effects can 
also be expected for p-wave pairing.ca' c^ Due to the boundary condition eq. ( 2.11 ) the gap equation 



(2.14) shows A2:=0 at x=0. This is related to the fact that the sign of the p^^-wave order parameter 
changes after the reflection at the surface. In Fig. 2 we show the self-consistently obtained order 



parameter, the current density, the magnetic field and the vector potential. In the bulk region A^, 
and Ay take the same magnitude as expected. Even if we start from an arbitrary relative phase, 
it finally becomes ib7r/2 yielding the pxUpy-waye state. In the GL theory this effect is included in 
the following fourth order term of the GL free energy which is derived in the Appendix, 



^4 2|A,|2|Aj,|2 + -(aX*' + A/'a2) . (3.4) 

The first term represents the competition between A^; and Ay, while the second term favors ±7r/2 
for their relative phase. Since the coefficient of the first term is four times larger than the second 
one, Px- and p^^-waves have to compete anyway. Near the boundary Ay is enhanced as shown in 
Fig. 2(a). This enhancement is connected with the suppression of A^^ near the surface, which opens 
the way for Ay to appear as a "single" order parameter. However, we find that A^ at the boundary 
(about 1.24) is larger than Asingic(O) (Asingic(0)/A(0)~1.21). Note that similar enhancement also 
appears in d-wave superconductors.E3'B3) The reason can be seen in the following sixth order term 
of the GL free energy, 

-4AQ\Ay\^\DxAx\^. (3.5) 

This term supports the A^-component whenever A^,. has a spatial variation and explains why 
Ay>Asingie near the surface. Microscopically the enhancement of Ay near the surface originates 
from the enhanced local density of states due to additional (bound) states appearing at low energy 
close to the surface as a result of the depletion of A^^ (as we will see below). Obviously only the 
Py-component can benefit from this additional density of states near the Fermi energy, since it does 
not suffer from pair breaking for the surface perpendicular to the x-axis. 

Due to the broken time-reversal symmetry, the local current fiows along the surface as shown 
in Fig. 2(b), whose direction is reversed if we change the relative phase of the order parameter 
components (px+'^Py'^-^Px — '^Py) (time reversal operation). The current along the x-direction is zero 
as expected. The magnetic field is induced by the surface current and it is screened by the Meissner 
effect, so that the total current becomes zero.Cj^ Note that depending on the geometry of the sample 
these currents can generate a finite magnetization. 

3.2 Density of states 

Due to the pair breaking effects we expect a modification of the local density of states near the 
boundary. The Andreev refiection and the conductance for the triplet pairing states were studied 
previously by various groups that the features of the conductance peak, which are related to the 
bound states, depend strongly on the angle of the incidence. E3'B3) The angle-resolved tunneling 
spectroscopy is related to the local density of states for a fixed Fermi momentum which is obtained 
from the quasi-classical Green function, 

N{ek,E,x) =Re g{9k,u;m,x) \iuj,r,-^E+i5 , (3.6) 



where (5 is a positive infinitesimal real number. Before examining the self-consistent order parame- 
ters case, it is illustrative to start with the consideration of a uniform order parameter, neglecting 
also the vector potential. In this case we can obtain the quasi-classical Green function analyti- 



g{Ok,^m,x) = gx{Bk-,^m,-,x)fx + gy{9k,UJm, x)fy + gz{6k,UJm,x)fz, 



9x{6k,UJm,x) = gx{TT - 9k,iOm,x) 

gy{Ok,^m,x) = -gy{TT - 6k,UJm,X 
gz{Ok,(^m,x) = gzilT - 6k,UJm,x) 



-A, 
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-A, 



1 + 



^m 



'2qx 



n 
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2 + A2( 

-2qx 



^m I — y V 



. ^m\^x{Qk)\ 



-2qx 



if^!A,(0fe)|A, 






y\^k) 2qx 



A^(0fc) = A(r) cos Ok, AyiOk) = A(T) sin 6^, n 



UJr, 



A2(r), q 



n 



VFx 



(3.7) 



Here fj {i=x, y, z) is the Pauli matrix in the charge space. We have assumed px+ipy-wa,ve symme- 
try, i.e. A{9k)=Ax{0k)+iAy{6k)=A(T){cos9k+ism9k), where A(T) is the magnitude of the order 
parameter in the bulk region at temperature T. At low temperatures we obtain A(T)~A(0). Note 
that the enhancement of Ay near the boundary in Fig. 2(a) can be understood by the first term 
of fe 



A, 



1 + 



Ai( 



m' + AK 



-2qx 



UJ. 



> A, 



(3.8) 



The vanishing of A^, at the boundary can be derived from the equation for gy . Substituting g^ into 
eq. ( |3.6| ), we have the following expression for the local density of states, 



N{ek,E,x)=N{TT-9k,E,x) 



rl-B| 
P 



A^i 



\E\ cos{2r]x) + sgn{E)\Axiek)\Ayiek)psm{27]x) 



+ 7r|A^(6'fc)|exp(^- 



s2_A2(r), 



2|A,( 



\VFx\ 

v = 



p[E^-Al{9k)] 

p 



A,, 



e(£;2_A2(r)) 



VFx 



(3.9) 



The first term of eq. (3.9) represents the continuum quasiparticle states, while the second term 
describes the bound state. For Aj^=0 the bound state is located at E=0 while a finite Py-wave 
part leads to E=Ay{9j^)=A{T) sm9kS^ This breaks the symmetry with respect to ky^—ky, so 
that the quasiparticle bound states with positive and negative ky have different occupation at low 
temperatures. Consequently, the quasiparticles generate a finite local surface current flows. For 
the uniform order parameter it is easy to see that Jy{x) has a finite value at x={) due to the jA^^jA^ 
term of ^^ in eq. (3.7), at T=0 Jj,(0)=efFA^(0)^A(0). This shows that the number of quasiparticles 
contributing to the current is of the order of A^(0)A(0). 

We now turn to the results of the full selfconsistent calculation. Figure 3 shows the local den- 
sity of states for a fixed momentum ky>{) {9k=iT/8) from the self-consistently determined order 



parameter and vector potential. The energy of the bound state estimated from eq. (3.9) is 
sin(7r/8)A(T)~0.38A(T). In the ky<0 case this energy is reversed in sign, since g has the sym- 
metry of eq. (2.7). Therefore the integrated density of states over the Fermi surface is symmetric 
under E^—E. In the case of the uniform order parameters the total local density of states at 
the boundary exactly becomes the same as the normal state, i.e. constant. One can test this by 
integrating eq. (3.9) over 9^ at x=0. We show the local density of states calculated from eq. (3.9) 
in Figs. 4(a) and (b). When we move away from the boundary, the bound states rapidly decreases 
and the bulk property is recovered. 

On the other hand, we notice gap structures in Fig. 4(d) even at x=0 for the self-consistent 
order parameter. At the surface we can see two gap energies, one is £^/A(T)~l, and the other 
forms a double peak feature around about 1.25. The first one corresponds to the bulk gap energy, 
while the second one is identified as the magnitude of Ay at x=0, which is enhanced near the 
boundary. These higher energy peaks are mainly formed by quasiparticles moving almost parallel 
to the surface as shown in Fig. 5(a), since they experience the enhanced Ay over a long distance (see 
Fig. 5(b)). According to eq. (0) the energy should be replaced by E—^E—evpyAy in the presence 
of a vector potential. This yields an energy shift in the density of states. Quasiparticles running 
almost parallel to the surface are strongly affected by the vector potential due to the factor vpyAy. 
Without the vector potential, the higher energy peaks merge into a single peak. Hence, the vector 
potential leads to the splitting into two peaks corresponding to the left and right moving particles 
along the surface. On the other hand, the bound states are mainly formed by the quasiparticle 
moving almost perpendicular to the surface shown in Figs. 5(c) and (d), which are only weakly 
affected by the vector potential. The self-consistently determined A^; shows a continuous spatial 
change near the surface which leads to a reduction of the magnitude of the delta function peak 
corresponding to the bound state in N{9k,E,x) compared to the uniform order parameter case. 
The reason is the widening of the potential well in which the bound state quasiparticle is trapped. 
The spatial dependence of A^; can be approximately expressed by tanh(x/(rfF/A(T))), where r 
controls the decay length of A^;. The magnitude of the bound state is then estimated bycJ) 

Kr) - ^(^ + ^/^) (3 10) 

^^'^^"r(i/2)r(i + r)- ^''•^^^ 

Here F is the gamma function, and I{r) takes 1 for r=0 and decreases with the increase of r. For 
the self-consistent order parameter r and I{r) can be estimated as 1 and 0.5, respectively. Therefore 
the density of states of the bound state at x=0 in Fig. 4(d) is about half of Fig. 4(b). 

Consequently, the local density of states near the boundary is sensitive to the spatial dependence 
of the order parameter, so that the self-consistent treatment is needed for quantitative discussions. 
However, the approximation of the uniform order parameters is helpful to capture the essence of 
the boundary effect. 



§4. External Magnetic Field 

In this section we study properties of the surface in an applied external magnetic field -Bext to 
the px+ipy-wave state. In this case the boundary condition for the magnetic field is Bz{0)=Bex_t, 
so that the first equation of (2.16) has to be replaced by 

B^ix) = 5ext - Ai r dx'Jyix'). (4.1) 

Jo 

We show the results of the self-consistent solution in Fig. 6. The order parameter does not suffer 
a significant modification compared with the situation without the external field. In Fig. 6(d) we 
can see that the peak position, which corresponds to the enhanced A^^, shifts in a higher (lower) 
energy region if we apply a negative (positive) field to the surface. The peak mainly comes from 
the quasiparticle having momentum of 6k~7r/2 as discussed before, which acquire an energy shift 
—evpyAy. In the case of negative B^xt, there is an additional positive contribution to the vector 
potential which shifts the peak related to the gap towards higher energies. If -Bgxt is reversed, 
the vector potential decreases and is almost zero at the surface when Bcxt=0-lBQ. Then the two 
peaks merge into a single peak as shown in Fig. 6(d). Increasing the external field further in the 
positive direction, the vector potential is reversed (see Fig. 6(c)), so that the peak splits again. The 
lines A and B in Fig. 6(d) correspond to the peak with the positive and negative ky, respectively. 
Energetically the case of the negative field is expected to be stable as it further opens the gap. 

In the d-wave case bound states are located at E=0f3 so that the peak structure at the zero 
energy is split by the vector potentialM^) In the present Px^ipy-^wave case bound states distribute 
inside the energy gap and form an almost flat local density of states at the surface as shown in 
Fig. 4(d). Near zero energy the shift is small due to the small vpy factor, since the states near zero 
energy are formed by the quasiparticles having momenta of 0^.^7r/2. Therefore the energy shift 
due to an external field is not observable for energies below the gap. The gap structure sensitive 
to the applied external magnetic field is located rather at energies corresponding to the enhanced 
Ay at the surface. 

§5. Domain Wall 

A modified quasiparticle spectrum is also expected near the domain wall. Clearly, domain walls 
cost energy and are not desirable modification of the order parameter. However, they exist even at 
low temperatures, since once they have formed at the onset of superconductivity, they are easily 
pinned at defects in the material and cannot move out of the sample. Domain walls between the 
pinning centers may give good targets for scanning tunneling spectroscopy on the sample surface. 
Physical properties of domain walls in unconventional superconductors have been investigated in 
the framework of the GL theory.SE^* Here we investigate their quasiparticle properties within the 
quasi-classical approximation which is essential for the tunneling spectroscopy. We will not examine 
effects of the external magnetic field in this case, because the magnetic field would not reach the 
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domain wall due to the Meissner screening. 

There are two kinds of domain walls separating the two degenerate superconducting states. 
For the domain wall perpendicular to the x-axis we have the {px—^Py\Px+'^Py) structure and the 
{—Px+^PylPx+Wy) structure, i.e. one of the two order parameter components should change sign at 
the domain wall. From the GL formulation it is expected that the first type is more stable, since 
the gradient free energy of the px-wave is three times larger than that of py-waye (see F^'^^ in eq. 
(A-11)). 

We assume that there is no scattering at the domain wall leading to quasi-classical trajectories as 
shown in Fig. 7. The boundary condition for the domain wall case is also given by eq. ( ^.11 ) with 



fcFi=^F2- III the domain wall case Jy{x) is symmetric under the x— >— x transformation. Therefore 
the vector potential is also symmetric and the magnetic field becomes anti-symmetric. 

5.1 The (px — iPy\Px+^Py) domain wall 

In Fig. 8 we show the results of the self-consistent solution. A current is induced near the domain 
wall and it is screened due to the Meissner effect .E3) Figure 8(a) shows that the relative phase 
between the px- and py-wave component suddenly changes, if we go through the domain wall (type 
I). On the other hand, the relative phase can change continuously, approximately as ^tanh(x/r,^o) 
(type II). These two types were discussed previously within the GL formulation£3E£P In the case 
of type II domain walls the relative phase changes gradually from +7r/2 to — 7r/2. Our calculation 
shows that the type I solution is energetically favored for the cylindrical symmetry Fermi surface 
case. 

Despite several similarities with the order parameter structure between the surface and the do- 
main wall, there are important differences in the quasiparticle spectrum as seen by comparing Figs. 
4(d) and 8(d). In Fig. 4(d) a gap structure can be seen at £'/A(T)~1.25, while there is no gap 
feature at the corresponding energy in Fig. 8(d). In Fig. 9 we show the schematic spatial depen- 
dence of the self-consistent order parameters, which a quasiparticle encounters along the classical 
trajectory. As we mentioned in §3, the higher energy peak for the surface case is mainly formed 
by the quasiparticles which run almost parallel to the surface. The same situation occurs also in 
the domain wall case, where the quasiparticle with momentum 6k~0 mainly contribute to the gap 
structure. Compared to the surface cases the distance in which the quasiparticle experiences the 
enhanced order parameter A^; is short (see Fig. 9(b)), so that the gap structure does not appear 
in the {px—^Py\Px+^Py) domain wall case. 

There is another characteristic point in Fig. 8(d), a small gap near the zero energy. The density 
of states near the zero energy mainly comes from the region where 0fc~7r/2 for the domain wall case 
(see Figs. 9(c) and (d)). The vector potential can generate a shift of the energy towards higher 
values. In the 0fc~7r/2 case the magnitude of the shift can be estimated as evp Ay {0)':^0.1, which is 
shown in Fig. 8(b). The small gap in Fig. 8(d) is about 0.1, which is consistent to this estimate. 
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In fact, neglecting the vector potential, we find that the gap structure in the small energy region 
has disappeared. 

A further difference from the surface case is the v-shaped density of states below the bulk energy 
gap. The origin of this feature is the gradual spatial dependence of Ay sin 0^ on the classical 
trajectory with 0fc~7r/2 (see Fig. 9(d)), which gives rise to additional bound states close to the 
bulk gap energy,EiP so that the density of states at low energies is shifted up to higher energy region. 

5.2 The (—Px+'Wy\Px+'Wy) domain wall 

In this case the system also favors the type I behavior for the relative phase. Here the situation 
is very similar to the surface case. (In the surface case A^^ changes its sign after the reflection at 
surface.) Therefore most of the properties are the same and differences between them appear only 
on the Fermi wave length scale. 

If the Fermi surface is not cylindrical, then the properties of the domain walls can be modified 
and a type II domain wall may be realized in some cases. S^ 

§6. Summary 

In this paper we analyzed the properties of the order parameter and the quasiparticles of a 
Pa;iiPj/-wave superconductor near the surface and at a domain wall. For this purpose we used 
a quasiclassical approximation which allows us to determine selfconsistently the current density, 
magnetic field and vector potential distribution besides the order parameter. 

The spatial dependence of the order parameter at the surface shows a suppression for the compo- 
nent p II n where n is the surface normal vector, i.e. the pa^-component is suppressed for the surface 
perpendicular to x-axis. The perpendicular component, on the other hand, is slightly enhanced. 
In the surface region we found bound states whose energy strongly depends on the position of the 
quasiparticle momentum on the Fermi surface. The average over the Fermi surface shows that this 
bound state yields a density of states which is essentially constant within the gap at the surface. 
As we move towards the bulk region these bound states gradually disappear and the complete bulk 
gap is recovered. We have also demonstrated that an external field (parallel to the c-axis) modifies 
the local density of states differently for the two degenerate states. Similar to the case of the time 
reversal symmetry breaking surface of the d-wave superconductor we find that the bound states 
carry a spontaneous current whose direction is opposite for the two degenerate superconducting 
bulk states. Furthermore an external magnetic field modifies the local density of states at the 
surface significantly only for energies slightly above the bulk energy gap, while the subgap state 
structure is essentially unaffected. 

The domain wall shows various similarities with the surface. One order parameter component is 
suppressed while the other is slightly enhanced. In contrast to the situation at surface, however, the 
suppressed and enhanced components are the ones with p _L n and p || n, respectively, which yields 
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the narrower (energetically more favorable) domain wall than in the opposite case. At the surface 
the choice is not available since the role of the order parameter components is entirely determined 
by the boundary condition due to the surface scattering. There are, in principle, two types of 
domain walls: In the type I domain wall the suppressed order parameter component vanishes in 
the center of the domain wall. For the type II domain wall this component introduces a phase twist 
which yields a finite modulus everywhere. The latter domain wall state is two-fold degenerate. The 
calculation using a cylindrical Fermi surface showed that the type I domain wall is more stable at 
all temperatures. This would change, if, for example, the anisotropy of the Fermi surface is taken 
into account. 

The local density of states in the center of the type I domain wall shows a tiny gap, while the 
structure of the large bulk gap remains visible. The shift of density of states is due to the presence 
of bound states in the domain wall similar to the surface and the tiny gap is a feedback effect of 
the vector potential created by the quasiparticles themselves. This change of the density of states 
could provide a possible way to detect domain walls by scanning tunneling microscopy. For a given 
voltage within the bulk gap the scanning of the (c-axis oriented) sample surface would yield an 
enhanced current in the region close the domain wall. 

Analogous to the surface case the bound states in the domain wall carry a spontaneous current 
which yields a magnetic field distribution. In contrast to the surface the net magnetization gener- 
ated by these currents at the domain wall has to be zero, since inside the superconductor magnetic 
flux can only enter in form of vortices enclosing a fixed quantum of flux. Thus, the magnetic fields 
of the domain walls are invisible on a macroscopic level. Even a scanning SQUID microscope has, 
at present, still too little spatial resolution to observe these fields. 

For the cylindrical Fermi surface the properties of the domain wall and the surface does not 
depend on the orientation of the normal vector as long as it lies within the x-y-plane. The real 
superconductor has certainly an anisotropic Fermi surface which could change various properties 
of the quasiparticle states for both the surface and the domain wall. To what degree changes occur 
remains a problem for future studies. 
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Appendix: Ginzburg-Landau Free Energy 

In §3 we have shown that the self-consistent order parameters favor the time-reversal symmetry 
breaking state (pxi^'ipy-wave) , and that the py-wave part is enhanced by the presence of the px- 
wave. In order to understand these results from a different point of view, we derive the GL free 
energy. The GL free energy can be easily derived by using the quasi-classical Green function.cJ' 
For simplicity the vector potential is neglected here. 

To derive the GL free energy, let us rewrite eq. (2.6) in the following form: 

U}mfikF,(^m,x) = A^kp, x)g{kp,U;m, x) - Df{kp,LL!m,x), 
U}mf{kF,(^m, x) = A*{kp, x)g{k-p,Ulm, X) + Df{kF,UJm, x), 

g^{kF,UJm,x) = 1 - f{kY,LOm,x)f{kY,UJm,x), (A-1) 

where D=^^-q^ and we have assumed the system is uniform in the y direction, i.e. there is no y 
dependence in 5, / and /. A(fcF,a;) is assumed to be composed oi Px- and p^^-waves. Since we are 
interested in the temperature near Tq, the order parameters are assumed to be small. The gradient 
term {D term) can be expressed as 

D _ \vyx\ d _^ \vyx\ d _^ d 

T ^ Tc dx^ A(0) dx ^ d{x/^o) ' ^ ' ^ 

where ^q=vf/ttA{0), which is the coherence length at r=0. For T^Tq, the spatial dependence of 
the order parameter is of the order of ^(T)'~t'F/vrA(T), and D /T^0{A/T)<^1. Let us treat A/T 
as a perturbation and derive the functions /, / and 5 in a power series of A/T, 

f = f0 + fl+f2 + ..., 

7 = 7o + 7i+72 + ---, 

9 = 90 + 91 + 92 + ■■■, (A-3) 

where the subscripts represent the order of A/T. The lowest order is contained in the bulk solution 

/bulk = / I lAI-J ' Jbulk = / ■) I I A I'i ' •'bulk = / '> , i a I'j ' (^'^j 

Therefore the zeroth order functions are given by 

/o = 0, 7o = 0, 90 = ^- (A-5) 

Since / and / are related to the order parameter, they do not have the zeroth order term. We 
can obtain the functions /, / and g up to any order of A/T by using eqs. (AT) and ( |A-5| )J1j The 
results are as follows: 

_ A 
^ A* 

= 0, 
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f2 
72 
92 ■ 
h 

53 - 



DA 
DA* 



lAp 



2aJm|tiJm| ' 

-A|Ap/2+D^A 

-A*|Ap/2+D^A* 

-ADA*+A*DA 

2aJm^|(^m| ■ 



(A-6) 



J. _ (-A^DA* + |ApDA)/2+DA|Ap/2-D3A 
j"4 — Til 3T7; i ' 

-J _ (A''^DA-|ApDA*)/2-DA*|Ap/2+D3A* 
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/s 



3|A|^/4-AD^A*-A*D^A+(DA)(DA*) 

2ulm^\uJm\ ' 



3|A|-*/4-AD2A*-A*D2A+{DA)(DA') 



-D 



-A2DA' + |ApDA+DA|A|2--2D3A 



2uJm'^\ljJm\ 

3|A|*/4-A*D2A-AD2A*+{DA*)(DA) -D 



-A*2dA+|A|2dA*+DA*|A|2-2D3A* 



2aJm'*|aJm| 



The gap equation can be written by using / and / as 

/A^.(x)' 

I 0<a;m<'^C 

Substituting / and / up to the fifth order into eq. ( |A-7| ), we obtain the gap equation for the 



A„fx)/ ' n.,^.,,„27r7-. l2sin(efc) 



f{9k,UJm,x)+ f {9kUJm,x) 



(A-7) 



Px-wave 






3A,IA,|2 + 2A,.|Aj,[2 + a:a2 



6A4l)'A, 



5A,|A,[4 + A^Af + 6A,|A,|2|Aj,|2 + SA^jA^fA^ + SA^lA^l^ + 2A*|Aj,|2A 



2 a2 

y 



2r.2, 



+5^6 A^L»^A; + 2A,IZ?,A,.1^ + 3A;(D,.A,)^ + 4|A,f D^,A,. 



lO^e^^A, 



+Ae[2AxAyDiA; + AAxAlDiAy + 2Ax\DxAy\' + AA^AyDiAy + SAHD^AyY + A^D^A^ 

+A\Ay\^DlAx + 2Ay{DxAx){DxA*y) + 2A^(D,A^)(D,.A*) + 6A;(D,.A,.)(£'.A^ 
, _, T 17 C(3) , _131 C(5) ^ _1 9 

and for the py-wave 



2 7~i2 A * 

0, 



-A2Ay + A4, 



3Aj^|Aj^P + 2Ay|A^P + A*A^ -2A4DlAy 



,3 A *2 



5A^|Aj,r + A;;Ar + 6Ay\Ay\'\Ax\' + 3A;|Aj^|"A^ + 3Aj,|A^r + 2A;|A^J"A, 



|2a2 



+Ae 



2Af,DtA 



x^y 



Af^DiA; + 2A,!Z?,A,|^ + SAliDxAyY + 4|A,|^Z)^A 
2AyAxDlAl + AAyAlDlA, + 2Ay\DxAx\^ + 4A;A,,Z?2a, + 3A;(Z?,A,)2 + AIdIa; 

= 0. 



+i\Ax\'DiAy + 2A,{D,Ay)iDxAl) + 2A,(Z),A,)(D,A:) + 6Al{D,Ay){D,A, 



(A-9) 



Here ^(2) is the Riemann zeta function. The GL free energy should recover the above gap equations 
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through the following Euler-Lagrange equation, 

d ^ d 



d 



dA* 



^"9(Z),A*)+^^a(I)2A*2) 



F = 0. 



(A-10) 



The result of the free energy density up to sixth order is as follows: 

F oc F2 + F4 + Fg + Fr'' + F^:^ + F^l^ + Fr'' 



6a 



6b 



6c ' 



F2 = -A2 

Fi = Ai 

Fa = -A 



|A,;r + |A,, 



^(lA.t^ + |A,|4) + 2|A.|2|A,p + ^(A^Af + Af A^) 

r: q Q 

A 1^ + lA l^') + -('IA ["^lA P + lA PlA 1^1 + -flA P + lA PVA^ A*^ + A*^A^) 



F^'^ = A, 



2 I ol n A |2 



Q\D^A.^\' + 2\D^Ay\ 



Fl 



rrad 



-A6i20\AJ'\a^AJ' + 



Al{D,Alf + Af{a,A 



1 



pgrad 
-^6b 



s;rad 



-A 



+A\Ay\'\D^Ay\' + - A'yiD^AlY + Al\D^Ay 



lQ\DlA^\^ + 2\D^Ay\^ 



6c 



-A 



1^ 

+2 

+4 



a2(d,a;)2 + Ai\D^Ayf + a2(z),a:)2 + a;^p,a, 
a,Aj,(d,a*)(d,a;) + a;a;(z),a,)(d,Aj,) 

A,A;(D,A,)P,A;) + AlAy{D^Al){D^Ay) 



+A,Al{D,Al){D,Ay) + A;Aj,(Z),A,)(D,A;) + \A,\^\D,Ay\^ + |Aj,|2[Z),A, 

(A-ll) 

The coefficient of the free energy density is chosen to reproduce the well known bulk GL free 
energy. The terms which suppress the spatial change of the order parameters are 6^4|D2:Aa;p and 
2A4\DxAy\'^ in the fourth gradient term. Therefore, the spatial variation of A^, is energetically 
more expensive than that of Aj^ along the x-direction. The asymmetry of the Px- and p^-waves 
comes from the D term, since D(xvfx=vf cos Ok- Fq^^ and F^^^ are the sixth order gradient free 
energies in which A^ does not couple with Ay, while Fq^^ is the free energy in which the two 
components are coupled. 
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Fig. 1. (a) Classical trajectory of a quasiparticle, in which the momentum of the incident and reflected quasiparticles 
along the surface is conserved, fcpi and feF2 are the momentum with feFx<0 and kFx>0, respectively, (b) ry is taken 
along the classical trajectory. A(fcFi,r||) and A(fcF2,r||) are the order parameters for the incident and reflected 
quasiparticles, respectively. A schematic spatial dependence is plotted. 
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Fig. 2. (a) Spatial dependence of the self-consistent order parameters for p^- and p^-waves. They are scaled by the 
bulk value A(T). The x coordinate is also scaled by ^o=t'F/7rA(0), where A(0) is the magnitude of the bulk order 
parameter at T=0. A^, and Ay are real and imaginary, respectively. The set of parameters are chosen as T=0.2Tc, 
uJc~WTc, k=Al/Co=2.5 (see ref. Eol). Here Xi^ — ^ym/e'^fin is the London penetration depth, (b) The self-consistent 
current density Jy{x), magnetic field Bz{x) and vector potential Ay{x). They are scaled by Jo=ewFA'^(0)rc, 
Bc='l'o/2\/27rCoAL and A(0)/euF, respectively. Here Be and ^o=h/2e are the critical magnetic field and flux 
quantum, respectively. 
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Fig. 3. (a) Local density of states for a fixed momentum, which has been calculated from the self-consistent order 
parameters. It is normalized by the one in the normal state. The energy E is scaled by A(T). The momentum is 
chosen to give Sfc=7r/8. A small imaginary part of 0.02 xA(r) is added to E for the plotting, (b) Local density of 
states at various positions. 
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Fig. 4. (a) Total local density of states for uniform order parameters. The vector potential is neglected for simplicity, 
(b) Total local density of states at various positions of (a) . (c) Total local density of states for the self-consistent 
order parameters, in which the vector potential is taken into account, (d) Total local density of states at various 
positions of (c). 
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Fig. 5. (a) A classical trajectory of an quasiparticle reflected by a surface, (b) A schematic spatial dependence of the 
self-consistent order parameters along the classical trajectory, ry is taken along the classical trajectory, (c) Same 
as (a) for a smaller 9k- (d) Same as (b) for a smaller 6k- 
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Fig. 6. (a) Spatial dependence of the self-consistent current Jy(x) which flows along the surface. Jy{x) and x are 
scaled by Jo=ei'FA''(0)rc and ^o, respectively. The parameters are chosen as T—0.2Tc, coc — WTc, k—2.5. (b) Spa- 
tial dependence of the self-consistent magnetic field Bz{x), which is scaled by the critical field Bc—^o/^^/^-k^o^l- 
(c) Spatial dependence of the self-consistent vector potential Ay{x), which is scaled by A(0)/e«F. (d) Local density 
of states at x=0 in an arbitrary unit. In each figure the applied external magnetic fields normalized by Be are 
depicted. In the preset case Bd is estimated as Bci/Bc=logK/\/2K~0.26. 
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Fig. 7. Classical trajectory of a quasiparticle going through a domain wall, in which the incident momentum is 
conserved. 
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Fig. 8. (a) Spatial dependence of the self-consistent order parameters for px- and pj,-waves. A^; and Ay are real and 
imaginary, respectively. The set of parameters are chosen as T=0.2rc, cjc^lOTc, k=2.5. (b) The self-consistent 
current density Jy{x), the magnetic field Bz{x) and the vector potential Ay{x). They are scaled by Jo=ei'pA'^(0)rc, 
Bc=^o/2V2-K^o^L and A(0)/e«F, respectively, (c) Total local density of states, (d) Total local density of states 
at various positions. They are symmetric under x^—x. 
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Fig. 9. (a) A classical trajectory of an quasiparticle going through a domain wall, (b) A schematic spatial depen- 
dence of the self-consistent order parameters along the classical trajectory, ry is taken along the classical trajectory, 
(c) Same as (a) for a larger 6^. (d) Same as (b) for a larger 9k- 
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